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Abstract 

We prove an approximation result showing how operators of the type 
-A - 7^(x - T) in L 2 (R 2 ), where T is a graph, can be modeled in 
the strong resolvent sense by point-interaction Hamiltonians with an 
appropriate arrangement of the 5 potentials. The result is illustrated 
on finding the spectral properties in cases when T is a ring or a star. 
Furthermore, we use this method to indicate that scattering on an 
infinite curve T which is locally close to a loop shape or has multiple 
bends may exhibit resonances due to quantum tunneling or repeated 
reflections. 



1 Introduction 

The main aim of this paper is to discuss a limiting relation between two 
classes of generalized Schrodinger operators in L 2 (K 2 ). One of them are 
measure-type perturbations of the Laplacian, in particular, we will be in- 
terested in potentials which are a negative multiple of the Dirac measure 
supported by a finite graph r C M 2 , in other words, they are formally given 
by the expression 

-A- 7 5(x-r) (1.1) 
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with some 7 > 0. We are going to show that such operators can be approxi- 
mated in the strong resolvent sense by families of point-interaction Hamilto- 
nians |AGHH| with two-dimensional S potentials suitably arranged. 

Stated in this way the problem is an interesting mathematical question; 
the indicated result will represent an extension the approximation theorem 
derived in jBFTj to the two-dimensional situation, which is not straightfor- 
ward because properties of the 5 potentials are dimension-dependent. In 
addition, we have a strong physical motivation coming from the fact that 
Hamiltonians of the type (jl.lj) are good models of various graph-type nanos- 
tructures which in distinction to the usual description [KSJ take quantum 
tunneling into account - see |Exlj for a bibliography to this problem. 

To investigate such models one has to find spectral properties of the op- 
erator (jl.lj) for various shapes of T. This is in general a complicated task 
even if the original PDE problems is refor mulated by means of the gener- 
alized Birman-Schwinger principle |BEKSj into solution of an appropriate 
integral equation. On the other hand, spectral problem for point-interaction 
Hamiltonians is reduced in a standard way to solution of an algebraic equa- 
tion with coefficients containing values of the free Green function. Hence an 
approximation of the mentioned type would be of practical importance; our 
second aim is to illustrate this aspect of the problem on examples. 

Let us describe briefly contents of the paper. In the next section we collect 
the needed preliminaries stated in a way suitable for the further argument. In 
Section|H]we formulate and prove our main result, which is the approximation 
indicated above. In the following two sections we discuss examples of two 
simple graph classes. The first are graphs of the form of a ring, full or 
open. In this case the spectrum of the operator (jl.lj) can be found by mode 
matching jETj which enables us to compare the approximation with the 
"exact" result, in particular, to assess the rate of its convergence. In contrast, 
in Section IS] we discuss star-shaped graphs. Here an alternative method for 
numerical solution of the spectral problem is missing, however, we can derive 
several conclusions for (an infinite) star analytically and show how does the 
approximation results conform with them. Computations of this type were 
performed already in jEJNj but the number of points used there was too small 
to provide a reasonable approximation. 

In the last section we address another question which is more difficult 
and no analytical results are available presently. It concerns the scattering 
problem on infinite leaky graphs with asymptotically straight "leads" . On a 
heuristic level, it is natural to expect that the states from the negative part 
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of the continuous spectrum can propagate being transversally confined to 
the graph edges, hence a nontrivial geometry should yield an Af x JV on-shell 
scattering matrix, where N is the number of the leads. This remains to be 
proved, however, and even more difficult it will be to compute the S-matrix 
mentioned above in terms of the graph geometry. 

For simplicity we will restrict ourselves to the simplest case when T is an 
infinite asymptotically straight curve, i.e. N = 2, without self-intersections. 
One natural conjecture is that if the distance between the curve points is 
small somewhere so that T is close to a loop shape locally, the system can 
exhibit resonances due to quantum tunneling. On the other hand, simple 
bends are unlikely to produce distinguished resonances; the reason is that 
the transverse 5 coupling in (jl.lj) has a single bound state, hence there is no 
analog here to higher thresholds which give rise to resonances in bent hard- 
wall tubes |DEMj . To support these conjectures, we have let us inspire by 
the so-called L 2 approach to resonances, which was in the one-dimensional 
case set on a rigorous ground in [HMJ. Specifically, we study the spectrum 
of a finite segment of T as a function of the cut-off position. In the first case 
we find in this dependence avoided eigenvalue crossings, the more narrow the 
smaller the gap to tunnel is, while no such effect is seen in the second case. 

On the other hand, a tunneling between arcwise distant points of the 
graphs is not the only source of resonances. To illustrate this claim, we 
analyze in our last example a Z-shaped graph with two sharp bends separated 
by a line segment. We find again avoided crossings the widths of which vary 
widely as functions of the bending angle. This provides an indirect but 
convincing indication that the above described on-shell S-matrix for a bent 
curve in the shape of a broken line is nontrivial. 

2 Preliminaries 

2.1 Generalized Schrodinger operators 

We start with the definition of Schrodinger operators of the form —A — 7m 
in L 2 (M. 2 ), where m is a finite positive measure on the Borel cr-algebra of 
T, which is assumed to be a non-empty closed subset of M 2 and also the 
support of the measure m. Furthermore, 7 is a bounded and continuous 
function, acting from V to K.+. We suppose that the measure m belongs 
to the generalized Kato class, in two-dimensional case it means that the 
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following condition holds, 

lim sup / | log(|x — y\)\m(dy) = , (2.1) 

e ^°x€R 2 Jb(x,e) 

where B(x,e) denotes the circle of radius e centered at x. 

Due to the fact that the measure m belongs to Kato class, for each a > 
there exists 6 G M such that any ijj from the Schwartz space 5(R 2 ) satisfies 
the inequality 

\ifj{x)\ 2 m{dx) < a \V^{x)\ 2 dx + b \^{x)\ 2 dx , (2.2) 
Jr 2 Jr 2 

as it was proven in the paper |SVj . Since 5(R 2 ) is dense in if 1 (IR 2 ) we can 
define a bounded linear transformation 



I m : H\R 2 ) h- L 2 (m) 

i m ip = V, W>gS(R 2 ). 

Using this transformation, the inequality (|2.2j) can be extended to the whole 
iJ 1 (]R 2 ) with the function ip on l.h.s. replaced by I m V ; - By employing the 
KLMN theorem, see |RS| Thm. X.17], we conclude that the quadratic form 
q im given by 

D{q im ) := H\R 2 ), (2.3) 
g 7m (^,0) := / Vtp(x)V(j)(x)dx - / I m -ip(x)I rn (j)(x)^(x)m(dx) , 



is lower semi-bounded and closed in L 2 (R 2 ). Hence there exists a unique 
self-adjoint operator H im associated with the form q im . 

The described definition applies to more general class of measures than 
we need here. If V is a graph consisting of a locally finite number of smooth 
edges meeting at nonzero angles, i.e. having no cusps, there is another way to 
define H im , namely via boundary conditions imposed on the wavefunctions. 
To this aim, consider first the operator H im acting as 



# 7m ) (x) = -{Ai/>)(x), xeR 2 \r 



for any ip of the domain consisting of functions which belong to H 2 (M. 2 \ T), 
are continuous at T with the normal derivatives having there a jump, 

P^{x) - p^{x) = ~^{x) , xer. (2.4) 
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Then it is straightforward to check that H im is e.s.a. and by Green's formula 
it reproduces the form q im on its core, so its closure may be identified with 
H im defined above. 

An important tool to analyze spectra of such operators is the generalized 
Birman-Schwinger method. If k 2 belongs to the resolvent set of H im we put 
i?* m := (Hyrn — k 2 )~ l . The free resolvent Rq is defined for Imk > as an 
integral operator with the kernel 

G k {x-y) = - A H^\k\x-y\). (2.5) 

Now we shall use Rq to define three other operators. For the sake of gen- 
erality, suppose that fi, v are positive Radon measures on R 2 with fi(x) = 
v(x) = for any x G R 2 . In our case they will be the measure m on V and 
the Lebesgue measure dx on R 2 in different combinations. By R* we denote 
the integral operator from L 2 (fi) to L 2 (y) with the kernel Gk, i.e. 

= °k * <i>t* 

holds z/-a.e. for all <p E D(R^) C L 2 (/i). 

With this notation one can express the resolvent R^ m as follows [B EKS : 



Theorem 2.1 Let Imk > 0. Suppose that I — jR^ nm is invertible and the 
operator 

from L 2 (R 2 ) to L 2 (IR 2 ) is everywhere defined. Then k 2 belongs to p(H im ) 
and {H irn - k 2 )~ l = R k . 

The invertibility hypothesis is satisfied for all sufficiently large negative 
k 2 because for such a k 2 the op erator norm of jR m!m (z) acting in L 2 {m) 



is less than 1, see again |BEKS , and similar result can be proven for the 
operator norm in L°°(m) following BFT . Thus from now on we consider 
k 2 < such that both these norms are less than one. 

For later considerations it is useful we rewrite operator H im = — A — 7m 
in the form —A — ^/i, where we have introduced 

/i = -j , a = -j . i^-oj 

J 7m J 7m 
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Since function 7 acquires only non-negative values and m is positive measure, 
a is a positive number. The resolvent reads 

/ 1 \ - 1 1 

It can be rewritten alternatively as 

(H, m - = + i&, M <7, (2.8) 

where t/> G L 2 (IR 2 ) and a G £ 2 (/i) represents the unique solution to the 
equation 

a* - = ^JU^ ^ - a.e. (2.9) 

Since -Rp^ belongs to if 2 (IR 2 ), using Sobolev's embedding theorem we con- 
clude that it has a bounded and continuous version, and therefore a also has 
a representative which is bounded and continuous on the set 5* 7 := {x G T : 
j{x) ^ 0}. 

2.2 Schrodinger operators with point interactions 

Consider a discrete and finite subset Y C T and the positive constant a 
defined above. As it is well known we cannot regard the operator Hy, a with 
the interaction supported by Y as before, i.e. as '—A + measure'. Instead, 
we define this operator via its domain: each function ip G D(Hy a ) behaves 
in the vicinity of a point a G Y as follows 

ip(x) = -— log \x - a\ L (ip, a) + L^ip, a) + 0(\x - a\), (2.10) 

Z7T 

where the generalized boundary values L (ip,a) and Li(^,a) satisfy 

Li(V»,a) + 2vr|F|aL (^,a) = (2.11) 

for any a G K with \Y\ := (1(^0; f° r a justification of this definition and 
further properties of the operators Hy i0l see, e.g., |AGHH| Chap. II. 4]. 

The form (j2.11|) which we have chosen is adapted for the indicated use of 
these operator in the approximation. In particular, the coupling parameter 
\Y\a depends linearly on the cardinality of the set Y . However, it is impor- 
tant to stress that our problem differs substantially form its three-dimensional 
analogue considered in jBFTl where no sign of a played a preferred role. In 
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the two-dimensional setting the coupling parameter tends to ±00 in the limit 
of weak and strong coupling, respectively. Consequently, we will be able to 
find an approximation for operators (jl.lj) with attractive interactions only. 

The Krein's formula for the resolvent (Hy, a — z)' 1 is the basic ingredient 
in the spectral analysis of the point-interaction Hamiltonians. 

Theorem 2.2 Let k 2 < and Ay jQ ,(/c 2 ) be the matrix \Y\ x \Y\ given by 



A Y , a (k 2 ;x,y) = — 
In 



2ix\Y\a 



Jxy 



-G k (x -y)(l- 8 xy ) 



(2.12) 



where Ce is the Euler constant. Suppose that this matrix is invertible. Then 
k 2 G p(Hy,o) and we have 



(H Y , a -k 2 )-^ (x) 



for any ip G L 2 



(x) 



[-^Y,a(k 2 )] ~ 1 (y, y') G k (x - y) 

y,y'& 

xR^iy') (2.13) 



One can see easily that the matrix Ay jQ ,(/c 2 ) is invertible for sufficiently large 



negative z 



-k . Indeed, the diagonal part is dominated by the term 



TjMogK, while the non-diagonal elements vanish as k — > 00, see the asymp- 
totic formula |AS[ 9.2.7] for the Hankel function -f^o 1 "*. In view of our special 
choice of the coupling, an alternative way how to make the matrix Ay )Q in- 
vertible is to take a sufficiently large set Y . 



Lemma 2.3 Let Imfc > and (Y T 

subsets of S 1 such that \Y n \ 
holds 

1 

sup — — sup 



n)neN be a sequence of non-empty finite 
00 as n — > 00 and the following inequality 



neN 



^2 G k(x -y) <a. 



(2.14) 

yeYn\{x} 

Then there exists a positive C and n G N such that the matrix Ay niQ ,(fc 2 ) is 

-1 



invertible and 



\Y, 



< c 



2.2 



(2.15) 



holds for all n > n . Here \\ ■ \\ p q means the norm of the map from £ p to 
the specification is superfluous here but it will be useful in the following. 
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Proof: Let us first decompose the matrix 1 / \Y n \ AY n , a (k 2 ) into the diagonal 
and non-diagonal parts, D n and R m respectively. For n being large enough 
the diagonal matrix D n is invertible and its operator norm in (C' y ™', || ■ H2) 
converges to a as n — > 00. Due to the strict inequality in the hypothesis there 
is an a < a such that the inequality (j2.14j) holds with a replaced by a. Then 
the Schur-Holmgren bound [AGHH| App. C] yields ||-R n ||2,2 < a < a, which 
in turn implies that the matrix sum D n + R n is invertible for a sufficiently 
large n. | 

In analogy with the expression ()2.8j) it is possible to rewrite Krein's for- 
mula f)2.13j) in the form 

(H Y , a - k 2 y l ip (x) = R k *P (x) + G *( x ~y)Qv, ( 2 - 16 ) 

y& 

where q y , y EY n solve the following system of equations, 

Qy~ E G k (y-y')q yl = ^(y), (2.17) 

for all y G Y n . 

3 Approximation by Schrodinger Operators 
with Point Interactions 

With the above preliminaries we can proceed to the main goal - we will prove 
that for a chosen generalized Schrodinger operator with an attractive interac- 
tion one can find an approximating sequence of point-potential Schrodinger 
operators under requirements which will be specified below. 

The assumption about positions of the point potentials is obvious - loosely 
speaking, as the sites of potentials are getting denser in the set T, they 
must copy the measure fi. Then we have to specify the coupling-constant 
behaviour in the approximating operators. We have already mentioned that 
in the analogous situation in dimension three the coupling constants scale by 
|BFTj linearly with the number \Y\ of point potentials, and suggested the 
same behaviour here. This requires an explanation, because it is well known 
that the coupling constants are manifested differently in dimension three and 
two; just consider a pair of point potentials and let their distance vary. 



1 

27 



2-|y|n-r]og(^)-r(\ 
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To get a hint that the scaling behavior for the approximation remains 
nevertheless the same, consider an infinite straight polymer as in [5QHH, 
III. 4], denote the coupling constant by a and the period by lon~ l . The 
threshold of the continuous spectrum is given as the unique solution E = — h 2 
to the implicit equation 

n 1 2im , v— > / 1 1 \ , s 
a = — log — — + lim > ; . 3.1 

21 k 2tt 6 l M ^^[ \ v / (27im) 2 + (K/o/n) 2 2nm J 

Now let the number n increase. If we want to keep the solution k preserved 
as n — > oo, then a must grow linearly with n; recall that for a > this means 
that the individual point interactions are becoming weaker. This motivates 
the choice of the coupling constants in the form \Y\a which we made in the 
boundary condition (|2.11j) . 

Now we can prove the announced approximation result. 

Theorem 3.1 Let V be a closed and non-empty subset of M? and let in be 

a finite positive measure on Borel a -algebra of T with supp m = T, which 
belongs to the Kato class. Let 7 : V —>■ M + be a nontrivial bounded continuous 
function. Choose k with Im k > such that the equation \2.9jl has a unique 
solution a which has a bounded and continuous version on S 1 := {x G T : 
7(x) 7^ 0}. Finally, suppose that there exists a sequence (Y n )^ =l of non-empty 
finite subsets of S~ such that \Y n \ — ► 00 and the following relations hold 



1 n| y&Y n J 

for any bounded continuous function f : T — > C, 
1 



su P'rrT su l ) Gk{x — y) < a , (3.3) 

V£Yn\{x} 



neN \Y n \ x &„ 



sup 



•* n.\ _ _ , , , 



yeY n \{x] 



(3.4) 



for n — > 00. The operators Hy n) a and H im defined in Section^ then satisfy 
the relation Hy n , a —* H jm in the strong resolvent sense as n — > 00 . 
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Proof: Since both Hamiltonians Hy n OL and H im are self-adjoint it is sufficient 
to prove the weak convergence, i.e. to check that 



?, {H Yn>a - z) 1 ip - (H jm - z) V)l2(r2) -> as n 



oo 



holds for arbitrary ip, <p G L 2 (M 2 ). Using the formulae f)2.8j) and f)2.16j) for 
the resolvents we arrive at 
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Y 



y'eY n 



The sum of the last two terms tends to zero as n — > oo, which follows from 
the hypothesis ()3.2|) . Since the function i?Q0 is continuous and bounded, as 
we have already mentioned, it is enough to prove the following claim, 

"f^lli —>■ as n — > oo 



for the £ x norm, where (u^)^ := 1^1?^ — a {y)i y £ ^n- 

To this end, we employ the equations (|2.9|) and (|2.17|) and we substitute 
the term R k dx ip in one equation from the other. The equation (|2.9|) holds \x- 
a.e., so we must consider continuous representatives of the functions involved 
here. In this way we get 



oc<y{y)-{Rl,»°){y) = ^ 



2vr«( Z /)|F n | + log( % ) +C E 



- Y G k{y - y')%' ■ (3-5) 

By adding two extra terms to both sides of the equation we arrive at 
1 Y-ta ,,2 ivi r '\m log(ifc/2) + C E 



Y [Ay n Ak 2 ;y,y')MY n \-a(y'))] = - ^;'^ B a{y) 
L — ' 27r|r n | 

\? Y G ^y - y>W - J G ^y - y>(y>(dy')- 



+ \Y 



We denote the vector on the r.h.s by w^ n \ then the previous formula reads 
(l/\Y n \)AY n ,a{z).v^ n > = w( n \ Since we assume that the inequality (jH.Hj) holds, 
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the Lemma 12.31 is applicable here. Therefore there exists no G N such that 
the matrix Tp-rAy^ jQ ,(z) is invertible for all n > n Q . Then we can write 



< 



\Y n 



\Y 



-1 



(n)| 



oo,l 



From the Lemma [2.31 and the relation UAH^i < \Y n \ \\ A\\ 2 ,2, A is a operator 
acting on C' y ™', we conclude that the operator norm in the inequality above is 
bounded by \Y n \C for some C > 0. Finally, let us look on the norm || H^: 
the first term of (w^ n ') y converges to zero uniformly w.r.t y as n — > 00 (recall 
that a is bounded) and the hypothesis ()3.4|) yields that remaining two terms 
go to zero as well. | 



4 Soft ring graphs 

Let us now pass to examples. The first class of graphs to which we apply 
the approximation developed in the Section El are rings, both full or open. 
Since the spectral properties of Hamiltonians with interaction supported by 
these ring graphs were already explored in the paper |ETj - see also the three- 
dimensional analogue discussed earlier in [ATjS] - this gives us an opportunity 
to compare the approximation with the "exact" results, in particular, to asses 
the convergence rate of the approximation. 

Consider a circle r := {x G M 2 : \x\ = R} with the radius R > and let 7 
be a function from r to R such that 7(2) = 7X[o,27r-e](v 9 )) where x = (R, (p), 
7 > and < 9 < 2ir. The Schrodinger operator with the 5 interaction 
supported by T given formally by (jl.lj) will be denoted by H^^; it can be 
given meaning in either of the two ways described in Sectional The spectrum 
and eigenf unctions are found easily in the full ring case, 8 = 0, when i? 7i i? 
is reduced by the angular momentum subspaces. Every eigenstate except 
the ground state is twice degenerate. In contrast, for a cut ring, 9 > 0, the 
spectrum is simple; the problem can be solved numerically by the mode- 
matching method |ET|. 

We start the presentation of the numerical results with the full ring. First, 
we put R = 10 and 7 = 0.5. The discrete spectrum of H 1>R consists of three 
eigenvalues, the ground state corresponds to the angular momentum I = 0, 
the other two correspond to / = ±1,±2, respectively, and they are twice 
degenerate; the number of eigenvalues is given be the inequality jR > 2\l\. 
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Figure 1: The dependence of eigenvalues of H Yol on the number of point 
potentials iV for 7 = 0.5 and the ring graph with R = 10. The dotted lines 
are the exact eigenvalues E = —0.0655, E\ = —0.0524, and E 2 = —0.0207. 

The choice of the approximating point-potential operators Hy n , a is obvious - 
N point potentials are spread periodically all over the circle and the coupling 
constant a equals N/(2ttRj). The task of finding the eigenvalues E of Hy n , a 
means to solve the implicit equation 



We plot the eigenvalues of H Y>a as iV = \Y\ increases in Fig. [TJ 

In case of a stronger interaction, 7=1, one gets a similar picture, just 
the number of levels rises to five, see Fig. El Note that the convergence of 
eigenvalues is slower than it is for the previous system. To estimate the rate 
of convergence, we calculate the difference between the exact eigenvalue and 
the eigenvalue computed using the approximation. The result is shown in 
Fig. |2J we observe that the aforementioned difference decays roughly as iV~ a 
with a being less than 1. 

The approximation by point-potential Schrodinger operators allows us to 
find easily the eigenf unctions. By |AGHHi Thm. II. 4. 2], they can be written 
as a linear combination of the free Green functions: 



detA Y>a ( E ) = 0- 



(4.1) 




(4.2) 
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100 200 300 400 500 600 700 800 900 1000 

N 

Figure 2: The dependence of eigenvalues of H Y , a on the number of point 
potentials N for 7 = 1 and the ring graph with R = 10 . The dotted 
lines are the exact eigenvalues E = —0.253, E\ = —0.243, E 2 = —0.21, 
E 3 = -0.159, and E 4 = -0.0881. 



4.5 5 5.5 6 6.5 7 

log(N) 

Figure 3: The dependence of the error on the number of point potentials N 
in the logarithmic scale. Dotted line corresponds to R = 10 and 7 = 0.5 and 
the solid line corresponds to R — 10 and 7 = 1. 

where k$ is the eigenvalue, i.e. detAy^/cg) = 0, and c is the solution to 
Ay, Q (fco) c = 0- The eigenfunctions obtained in this way behave as one ex- 
pects: they decrease exponentially if moving transversally away of the circle, 
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-2 8 

Figure 4: A detail of the wavefuction near the intersection of the circle and 
one of the nodal lines. It is the wavefunction of the third excited state (/ = 3) 
for R = 10, 7 = 5 and \Y\ = 100. 




and for I > they copy a sine function if moving along the circle. A closer 
inspection of the eigenfunctions (|4.2|) shows, of course, a logarithmic peak at 
the site of each point potential, as pictures Fig. 0] and Fig. demonstrate. 
To our opinion the contributions to energy coming from these spikes are 
responsible for the slow convergence of the approximation. 
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Figure 6: The wavefunction of the fifth excited state of H l R for R = 10, 
7 = 1 and 9 = n/3. The solid line represents the ring T. 




Figure 7: The wavefunction of the fifth excited state of Hy, a for \Y\ = 1000, 
which approximates the fifth excited state of H lt R from Fig. El The solid 
lines represent the nodal lines, the dotted line represents the graph V. 

Next, we consider an open ring, 9 = tt/3, i.e. one sixth of the perimeter 
is missing. For example, the fifth excited state for R = 10 and 7 = 1 has 
the energy E 5 = —0.151, the corresponding eigenfunction is shown in Fig. El 
The approximation by 1000 point potentials yields the energy E' 5 = —0.116 
and the corresponding eigenfunction shown in Fig. [71 
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5 Star-shaped graphs 



Another class of leaky-graph systems to which we will apply the approxi- 
mation by point-interaction Hamiltonians are the star-graph Hamiltonians 
where T is a collection of segments coupled at a point. In distinction to the 
previous section no "direct" method to solve the problem is available in this 
case, and thus the approximation represents the only way how to obtain a 
numerical description of the eigenvalues and eigenfunctions. 

First we want to make a remark about fmiteness of such graphs. Our 
interest concerns primarily infinite stars in which the arms are halflines, in 
particular, since they support localized states despite the fact that the graph 
geometry would allow escape to infinity. On the other hand, the above ap- 
proximation result applies to finite graphs only because it requires J 7m to 
be finite. Nevertheless, the result can be used, due to the fact that the infi- 
nite star Hamiltonian is approximated, again in the strong resolvent sense, 
by a family of operators with cut-off stars. To justify this claim, it is suffi- 
cient to realize that the corresponding family of quadratic forms is by (|2.3|) 
monotonous and bounded from below, so that Thm. VIII. 3. 11 of |Kaj applies. 

Let us thus discuss in the beginning what can be derived analytically 
about infinite star graphs. Given an integer N > 2, consider an (N— l)-tuple 
(3 = {/3i, . . . , Pn-i} of positive numbers such that 

iV-l 

(3 N := 2tt- > °- 

Denote $j := 5^i=i/%> where we put conventionally i?q = 0. Let Lj be the 
radial halfline, Lj := {x G 1R 2 : argx = $j}, which is naturally parametrized 
by its arc length s = \x\. The support of the interaction is given by T = 

r(/?):=Uf= % 

The star-graph Hamiltonians H^{I3) are defined formally by (jl.ljl as 
Schrodinger operators with an attractive potential supported by the graph 
T and with the coupling constant 7, and the proper meaning is given to this 
operator in the way described in Section |21 

Remark 5.1 Properties of the operator H^{f3) certainly depend on the or- 
der of the angles in (3. However, the operators obtained from each other 
by a cyclic permutation are obviously unitarily equivalent by an appropriate 
rotation of the plane. 
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Let us mention two trivial cases: 

Example 5.2 (a) ^(tt) corresponding to a straight line can be written 
as h y ® I + I ® (— <9y), where /i 7 = — <9 2 — j8(x) is the one-center point- 
interaction Hamiltonian AGHHJ on L 2 (R). Consequently, its spectrum is 
purely a.c. and equal to [— 7 2 /4, oo). 

(b) H±((3 S ) with (3 S = {f,f,f} again allows separation of variables being 
/i 7 <E> I + I <8> h 7 . Hence the a.c. part of a(H^(P s )) is the same as above, and 
in addition, there is a single isolated eigenvalue — 7 2 /2 corresponding to the 
eigenfunction (2 7 )- 1 e - 7(|a;|+ls/|)/2 . 

5.1 The essential spectrum 

First we notice that the essential spectrum of H^{(5) does not extend below 
that of -f^ 71 ") corresponding to a straight line. 

2 

Proposition 5.3 inf 0" e ss(-^iv(/?)) > — \ holds for any N and j3. 

Proof: By Neumann bracketing. We decompose the plane into a finite union 



where Wj is a wedge of angle (3j, Sj is a halfstrip centered at Lj which is 
obtained by a Euclidean transformation of M + x [£, £] , and P is the remaining 
polygon containing the vertex of T. Imposing Neumann boundary conditions 
at the common boundaries, we get a lower bound to H N (f3) by an operator 
which is a direct sum corresponding to the decomposition 1)5.1)1 . The wedge 
parts have an a.c. spectrum in IR + while the polygon has a purely discrete 
spectrum. Finally, the halfstrip spectrum is a.c. again and consists of the 
interval [e(£), oo), where e(£) is the lowest eigenvalues of (— <9 2 — j5(y))N on 
L 2 {[— £,£]). It is straightforward to see that to any 77 < — 7 2 /4 there is £ 
such that e(£) > 77, and since the decomposition (j5.1)l can be chosen with £ 
arbitrarily large, the proof is finished. 1 

In fact, the essential spectrum is exactly that of a straight line. 

Proposition 5.4 ct css (Hn((3)) = [— 7 2 /4, 00) holds for any N and (5. 




(5.1) 
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Proof: In view of the previous proposition, it is sufficient to check that 
a ess (H N (f3)) D [-7 2 /4,oo). Given a function e C °°([0,oo)) with ||0|| = 1 
and </>(r) = 1 in the vicinity of r = 0, we construct 

ip n (x;p,x n ) := —^==<P (-\x-x n \) e -^ (2) \/^ e ^ m 

with p > and x = (x^\ a^ 2 -*) , where the points be chosen, e.g., as 

(n 2 , 0). It is easy to see that the vectors if) n — > weakly as n — >■ oo and that 
they form a Weyl sequence of Hn(/3) referring to the value — 7 2 /4 + p 2 . This 
yields the desired result. 1 



5.2 The discrete spectrum 

The first question concerns naturally the existence of isolated eigenvalues. It 
follows from two observations of which one is rather simple. 

Proposition 5.5 H^{f3) > Hn+i{(3) holds for any N and angle sequence 

$ = {/?!, • • • , /%-!, Pf\pf\fr +li /3v-l} With + ft? = ft . 

Proof follows directly from the definition by the quadratic form (|2.3jl . I 

On the other hand, the second one is rather nontrivial. It has been proven 
in jEl] for a wide class of piecewise continuous non-straight curves which 
includes, in particular, a broken line. 

Proposition 5.6 CdiscO^O^)) ^ s nonempty unless (3 = it. 

Combining these two results with the minimax principle (recall that H^(f3) 
is below bounded) we arrive at the following conclusion. 

Theorem 5.7 cr disc (i?jv(/3)) is nonempty except if N = 2 and (3 = tt. 

Next one has to ask how many bound states does a star graph support. 
The answer depends on its geometry. There are situations, however, where 
their number can be large. 

Theorem 5.8 Fix N and a positive integer n. If at least one of the angles 
f3j is small enough, card (<7disc {H n (/?))) > n. 
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Proof: In view of Proposition 15.51 it is again sufficient to check the claim for 
the operator H 2 (P). Choose the coordinate system in such a way that the 
two "arms" correspond to arg# = ±/3/2. We employ the following family of 
trial functions 

*(x,y) = f(x)g(y) (5.2) 

supported in the strip L < x < 2L, with / £ C 2 satisfying f(L) = f(2L) = 0, 
and 

1 ... \y\ < 2d 

e -j(\ y \-2d) _ _ _ \y\>2d 

with d := Ltan(/3/2). Let us ask under which conditions the value of the 
shifted energy form 



g{y) 



*]Hiv*ii 2 -^ii/r4w 

is negative. Since ||g|| 2 = 4d + 7 _1 and \\g'\\ 2 = 7, this is equivalent to 

2 8ecf- 7 d-f 



1 + 4^d 

By minimax principle the system there will be at least n isolated eigenvalues 
provided 

7m /A 2 . \\fT ^ 2 2secf- 7 rf-| 

tan — = inf sup < 7 



d 2 J M ,i M„ \\f\\ 2 ' l + ^d ' 
where M n means an n-dimensional subspace in L 2 ([L,2L]), i.e., if 



TT 2 V 2 / 1 + 4 7 d 

Now one should optimize the r.h.s. w.r.t. 7c?, but for a rough estimate it is 
sufficient to take a particular value, say 7c? = sec f — | which yields 



1 (3 (8secf-5) 3/2 

— cot - - '- 

167T 2 ( 8sec |-3) ] 



n < cot — - 2 ; (5.3) 



it is obvious that the last inequality is for any fixed n satisfied if (3 is chosen 
small enough. 1 
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Corollary 5.9 Independently on (3, card (crd isc (H N (P))) exceeds any fixed 
integer for N large enough. 

Remark 5.10 The estimate used in the proof also shows that the number 
of the bound states for a sharply broken line is roughly proportional to the 
inverse angle, 

o3/2 

as (3 — > 0. This is the expected result, since the number is given by the length 
of the effective potential well which exists in the region where the two lines 
are so close that they roughly double the depth of the transverse well. 

5.3 The Birman-Schwinger approach 

Now we are going to show how the spectral problem for the operators H^{j3) 
can be reformulated in terms of suitable integral operators. We will employ 
the resolvent formula for measure perturbations of the Laplacian derived in 



BEKS . Notice that this technique was used in [El] to derive a result which 
implies our Proposition 15.61 

Since the operators H^{(3) are defined by the quadratic form ()2.3|) they 
satisfy the generalized Birman-Schwinger principle. If k 2 belongs to the re- 
solvent set of HjfiP) we put R* r := (H N (j3) — k 2 )' 1 . We already know the 
Krein-like formula for the resolvent from Thm 12.11 here it has the form 

with m denoting the D irac me asure on T. One can express the generalized 
BS principle as follows [BEKS : 

Proposition 5.11 dimker (H^(f3) — k 2 ) = dimker(J — 7_R^ m ) for any k 
with Im k > . 

Consequently, the original spectral problem is in this way equivalent to 
finding solutions to the equation 

K 7 = ^ ( 5 ' 4 ) 

in L 2 (T), where TZ* T := . Furthermore, in analogy with A.GHH, 

Sec. II. 1] one expects that a non- normalized one corresponding to a solution 
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4> of the above equation is 

ip(x) = J G iK (x-x(s))(j)(s) ds ; 
a proof of this claim can be found in the sequel to the paper |Poj . 

5.4 Application to star graphs 

Let us look now how the equation (J5.4j) looks like for graphs of the particular 
form considered here. Define 

dij(s, s') = dg(s, s') = ^s 2 + s' 2 -2ss'cos\$ j -^| (5.5) 

with -&j - -&i — Yji=i+\l3h in particular, du{s,s') — \s - s'\. By 7^(/3) = 
HjiiP) we denote the operator L 2 (IR + ) — » L 2 (R + ) with the kernel 

then f)5.4j) is equivalent to the matrix integral-operator equation 

TV 

£(W)-V)& = 0> i = l,...,N, (5.6) 

i=i 

on ®f = i L 2 (M. + ). Notice that the above kernel has a monotonicity property, 

7^(/3)>7^.(/3') (5.7) 
if — < — This has the following easy consequence: 

Proposition 5.12 Each isolated eigenvalue e n (/3) of H 2 (/3) is an increasing 
function of (3 in (0,7r). 

Proof: Notice first that the eigenfunction related to e n (/3) is even with respect 
to the interchange of the two halflines. Without loss of generality we may 
assume that argL,- = (— 1)- ? ~ 1 /3/2. The odd part of H 2 (f3) then corresponds 
to the Dirichlet condition at x — 0. The spectrum of this operator remains 
the same if we change argLi to 7r — (3/2. Removing then the Dirichlet 
condition, we get the operator H 2 (tt) with inf a{H 2 {^)) = inf o- css (H 2 (j3)) , so 
by minimax principle no eigenfunction of H 2 (/3) can be odd. 
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On the symmetric subspace the diagonal matrix element of the matrix 
integral operator TZ K (/3) equals 

(0,^03)0) = 2(0^00 +2(0!, K«M<l>i) 

with the first term at the r.h.s. independent of 3. Next we we notice that 
{7Z K ((3)} is a type (A) analytic family around any (3 G (0,7r); the derivative 

-^1Z K (f3) is a bounded operator or the form ^ ^ ^ ^ where Vp has the 

kernel 

^ / ,n 7 T , / ,/? / /n\ ss' sin 3 
Vp(s, sf) = -l K, Udf 2 (s, s')) jn-z- < o . 

"12 I s ? S ) 

At the same time {1Z K (f3)} is a type (A) analytic family w.r.t. k around 
any k G (0, oo) and the corresponding derivative is a bounded operator 

( S' 1 S' 12 ) with the kernel 

V 'M2 '^22 / 

7^(8, = " ^ [*o («dg(a, *')) + «dg(a, O^i *'))] < . 

Let 0^ = ( I) be a normalized eigenvector of lZ K (j3) corresponding to the 
eigenvalue \(k,@) = (cj) 13 ) TZ K (f3)4> 13 ). By Feynman-Hellmann theorem we find 



d 

d/3 



A(«,/3) = 2(0f,P^) <0, 



and similarly j- X(k,(3) < 0. The solution /t = of the implicit equation 
X(k, 3) = 1 thus satisfies 

and consequently, the eigenvalue — k((3) 2 is increasing w.r.t. (5. 1 



5.5 Numerical results 

Having explained analytically how does the discrete spectrum of H N (f3) de- 
pend on the number of arms and the angles /3, we employ now the approxima- 
tion by point-potential Schrodinger operators to obtain the numerical results 
for cut-off stars which would illustrate the above conclusions. 
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Figure 8: The dependence of the discrete spectrum of Hy >a on the angle (5 
for the symmetric two-arms star T\. The dotted line represents the threshold 
-7 2 /4- 

Consider a two-arms star graph r\ with both arms having the same 
length, L\ = L2 = 300, and put 7 = 0.1. As we already know, the threshold 
of the continuous spectrum of H 2 (P) is — 7 2 /4 = —0.0025. We approximate 
H 2 {f3) by point-potential Schrodinger operator Hy a that has 1 potential on 
the center of the star graph and 200 equidistant point potentials on each arm. 
Only the lower part of the discrete spectrum of Hy a approximates the dis- 
crete spectrum of the star-graph operator, while the upper part corresponds 
to the interval [— 7 2 /4,0] C a css (H 2 (/3)) . Therefore only the states with en- 
ergy below the threshold may be indeed understood as the approximation of 
the bound states of H 2 {(3). 

The main result of the analytic argument presented above was the de- 
pendence of the eigenvalues on the angle (3: if decreases, the eigenvalues 
decrease and their number grows. The discrete spectrum of Hy, a for (3 vary- 
ing is in a good agreement with this fact as Fig. |H] illustrates. The eigenvalue 
crossings we observe are actual, which is a consequence of the symmetry of 
the graph T\. For a non-symmetric graph the crossings become avoided. To 
demonstrate it in Fig. |UJ we slightly change the length of one arm, L 2 = 306, 
while the other parameters are preserved. Also the setting of the point po- 
tentials approximating the new graph T 2 stays the same, up to extra four 
potentials added on the longer arm. 

The expression (J4.2j) again yields the eigenf unctions in the form of a sum 
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of the free Green functions. We limit ourselves to a pair of examples: the 
ground state of H 6 (j3) in Fig. [TU] and the third excited state of H W (P) in 
Fig. El The length of the cut arms is 30, in the former case the approximat- 
ing operator Hy, a has 601 point potentials and in the latter case the number 
of point potentials is 1001. We see, in particular, that for N large enough 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 



Figure 9: The dependence of the discrete spectrum of H Ya on the angle 
/3 for the non-symmetric two-arms star T 2 . The dotted line represents the 
threshold — 7 2 /4. 




Figure 10: The wavefunction of the ground state of H Ym E = —0.612, 
which approximates the ground state of H 6 (j3) with (3 being the 5-tuple of 
7r/3 and 7 = 1. The contours correspond to logarithmically scaled horizontal 
cuts, the dotted lines represents the graph V. 
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Figure 11: The wavefunction of the third exited state of Hy^-, £3 = —0.265. 
It approximates the third excited state of H 10 ([3) with (3 being the 9-tuple 
of 7r/5 and 7 = 1. The convention is the same as above, the solid line is the 
nodal line. 

Hn(P) may have closed nodal lines; it is an interesting question what is the 
minimal N for which this happens. 

6 L 2 -approach to resonances 

Our last example, as indicated in the introduction, concerns the situation 
when T is a single infinite curve; we want to see whether resonances in the 
scattering of states propagating along T may be detected by inspecting the 
spectrum of the cut-off problem with the curve of a finite length which is a 
parameter to be varied. 




Figure 12: An example of the curve T: R = 10 and A = 1.9. 
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Figure 13: The dependence of the discrete spectrum of the Hy, a on the length 
L for R = 10, 7 = 1 and A = 1.9. 
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Figure 14: The dependance of the discrete spectrum of Hy a on the length L 
for R = 10, 7 = 1 and A = 2.9. 

Consider a curve T from Fig.^J its central part consists of three segments 
of a circle with the radius R = 10 and it has two infinite "legs". The distance 
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Figure 15: The dependence of the discrete spectrum of Hy, a on the length L 
for R = 10, 7 = 1 and A = 5.2. 

between two closest points of the curve, i.e. the bottleneck width, is denoted 
by A. The coupling constant 7 equals 1. We cut the "legs" of the graph to a 
finite length L and we plot the eigenvalues computed using the approximation 
for L varying. The number of the point potentials involved in Hy a is chosen 
so that the distance between every two adjacent points equals 0.3. For small 
values A = 1.9 and A = 2.9 the tunneling effect occurs and we can see the 
plateaux in Figs. [Till EH which indicate the existence of resonances; the width 
of the avoided crossings increase with A as expected. On the other hand, for 
more open curve, A = 5.2, we get a different picture, where the 'plateaux are 
absent, see Fig. HH1 

The second graph type we are interested in are simple bends, in terms 
of Section El it is a two-arms star T(j3). As we have already mentioned, the 
avoided eigenvalue crossings are not expected here because the transport 
along the graph arm involves a single transverse mode - this is confirmed in 
Fig- EE which shows the results of the cut-off method for f3 = 7r/4. 

Finally to illustrate that resonances may also result from multiple reflec- 
tions rather than from a tunneling, we apply the cut-off method to graphs 
which are of a stair-type, or Z-shaped. They consists of three line segments: 
the central one has a finite length R, i.e. the "heigth" of the stair, the other 
two are parallel (cut-off) halflines. If the angle 9 between the segments is less 
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Figure 16: The dependence of the discrete spectrum of Hy, a on the arm 
length L for the angle (3 = 7r/4 and 7 = 1. 

than 7r/2 the term Z-shaped is more appropriate. The results for R = 10 and 
7 = 5 are plotted in Figs. El and in the former case the stair is "skewed" 
to 9 = 0.327T, in the latter we have 9 = n; the distance between the point 
potentials of the approximating operator Hy, a equals 0.1. We find avoided 
crossings, however, they are very narrow for 9 = 0.327T and barely visible in 
the right-angle case; this observation can be naturally understood in terms 
of the reflection probability through a single bend - compare with the angle 
dependence of the spectrum on Fig. |H1 
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Figure 17: The dependence of the discrete spectrum of Hy t a on the arm 
length L for R = 10, 9 = 0.32vr and 7 = 5. 
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